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Abstract 

We include the full second-order corrections to the static QCD potential in the analysis 
of the tt threshold cross section. Then we examine the difference between the results 
obtained in the momentum-space approach and in the coordinate-space approach, which 
was found recently. Contrary to our expectation, the reduction of the difference by the 
inclusion of the second-order corrections is very small. There still remains a non-negligible 
deviation, which originates from the difference in the construction of the potentials in the 
two spaces. We scrutinize this problem. In particular, we estimate our present theoretical 
uncertainty of the tt threshold cross section at the peak to be foipeak/cpeak ~ 6% within 
perturbative QCD. 



*On leave of absence from Department of Physics, Tohoku University, Sendai 980-77, Japan. 



In this paper we report on our present theoretical understanding of the tt total cross sec- 
tion near the threshold. Up to now, all the 0(a s ) corrections (also leading logarithms) have 
been included in the calculations of various cross sections near threshold. In order to take 
into account the QCD binding effects properly in the cross sections, we have to systematically 
rearrange the perturbative expansion near threshold. Namely, we first resum all the leading 
Coulomb singularities ~ (a s //3) n , take the result as the leading order contribution, and then 
calculate higher order corrections, which are essentially resummations of the terms ~ a™ +1 //3 n , 
a™ +2 / f3 n , ... It is also important to resum large logarithms arising from the large scale difference 
involved in the calculation JT|.[] This is achieved by (first) calculating the Green function of the 
non-relativistic Schrodinger equation with the QCD potential 0, []J. Conventionally both the 
coordinate-space approach developed in Refs. [|T], |3] and the momentum-space approach devel- 
oped in Refs. ||, f| have been used in solving the equation by different groups independently. 
It has recently been found P that there are discrepancies in the results obtained from the two 
approaches reflecting the difference in the construction of the potentials in both spaces. It was 
argued that the differences are formally of 0(a 2 s ) but their size turns out to be non- negligible. 

Quite recently there has been considerable progress in the theoretical calculations of the 
second-order corrections to the cross section at threshold and the Coulombic bound-state prob- 
lem. New contributions have been calculated analytically |7], [8] and numerically for QED 
bound-states. Very important steps have been accomplished in QCD as well. The full second- 
order correction to the static QCD potential was computed in f|3T0| | . Also, the 0(a 2 ) total cross 



section is known now in the region « s < /? < 1 as a series expansion in {3 |pH| . All these 
results have to be included in the calculation of the full 0(a 2 ) corrections to the threshold 
cross section, which has been completed just in these days (as far as the production process of 
top quarks are concerned) . 

In this report, we incorporate the full 0(al) corrections (the second-order corrections to 
the leading contribution) to the static QCD potential into our analyses. In principle this is 
a step towards an improvement of the theoretical precision in our analysis of the tt threshold 
cross section. The reduction of the difference is very small, however, and there still remains a 
non-negligible difference. We scrutinize this problem of the difference between the momentum- 
space and the coordinate-space potentials. We find that there is a theoretical uncertainty 
within perturbative QCD which limits our present-day theoretical accuracy of the threshold 
cross section. 

Let us first state the numerical accuracies attained throughout our analyses. We confirmed 
that our numerical accuracies are at the level of 10~ 4 . We have tested our programs with the 
Coulomb potential whose analytical form is known both in momentum space and in coordinate 
space. Moreover we confirmed that we obtain the same cross section within the above accuracy, 
irrespective of whether we solve the Schrodinger equation in momentum space or first Fourier 
transform the potential and solve the Schrodinger equation in coordinate space. In this way 
we also checked that our numerical Fourier transformation of the potential (from momentum 
space to coordinate space) works within the above quoted accuracy. The level of accuracies is 
quite safe in studying the size of the higher order corrections which are described in this paper. 



* Since the toponium resonance wave functions have wide distributions ~ 10-20 GeV, they probe a fairly 
wide range of the QCD potential. For example, this is reflected in the fact that the fixed-order calculation with 
any single choice of scale /i cannot reproduce simultaneously both the distribution and the normalization of the 
differential cross section which includes all the leading logarithms. It is known that the normalization of the 
cross section is more sensitive to the short-distance behavior of the QCD potential. 
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Let us now briefly explain the construction of our potentials in momentum space and in co- 
ordinate space, respectively. More detailed descriptions including formulas are given in the ap- 
pendices. The large- momentum part of the momentum-space potential Vjkpt(<z) is determined 
as follows. First the potential has been calculated up to 0(ofy in a fixed-order calculation. 
The result is then improved using the three-loop renormalization group equation in momentum 
space. At low momentum, the potential is continued smoothly to a Richardson-like potential. 
On the other hand, the short-distance part of the coordinate-space potential Vsfhmn ( r ) is cal- 
culated by taking the Fourier transform of the fixed-order potential in momentum space, and 
then is improved using the three-loop renormalization group equation in coordinate space. At 
long distance, the potential is continued smoothly to a phenomenological ansatz. Thus, it is 
important to note that the two potentials are not the Fourier transforms of each other even 
in the large-momentum or short- distance region. They agree only up to the next-to-next-to- 
leading logarithmic terms of the series expansion in a fixed MS coupling. The difference begins 
with the non-logarithmic term in the three-loop fixed-order correction. 
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Figure 1: Comparison of the total cross sections (normalized to R) calculated from the different 
potentials: Vjkpt (solid), Vsfhmn (dashed), and V new (dotted line). We set a^g(M|) = 0.115, 
m t = 175 GeV, and r t = 1.427 GeV. 

In Fig. |I] we show a comparison of the total cross sections (normalized to R) calculated from 
Vjkpt (solid) and from Vsfhmn (dashed line), without any weak or hard-gluon corrections: 

4 roo 

R=— dpp 2 \G(E,p)\ 2 T t . (1) 
nrnf Jo 

For the physical parameters we used a^jg(M§) = 0.115, m t = 175 GeV, and T t = 1.427 GeV. 
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We find that the two cross sections differ by 6.7% at the peaks and by 1.9% at E = 5 GeV.jj] 
Since the difference of the cross sections calculated from the next-to-leading order potentials is 
8.2% at the peak and 2.2% at E = 5 GeV for the same value of a s (M|) (see Figs. ||and|6|), the 
cross sections have come closer only slightly after the inclusion of the second-order correction 
to the potential. The remaining difference is much larger than what one would expect from an 
0(ot Z g) correction relative to the leading order, which is not fully included in our analyses, even 
if we take into account the high sensitivity to the coupling, cr pea k oc a 2 [|| . The purpose of this 
report is to understand the origin of this unexpectedly large difference. 

As already mentioned, the difference of the cross sections reflects the difference of the poten- 
tials. The derivative of the potential dV(r)/dr is directly related to the size of the cross section; 
the cross section is larger if dV(r)/dr (= magnitude of the attractive force) is larger. This is 
because, with increasing probability that t and t stay close to each other, the wave function at 
the origin (^(O)] 2 increases, and so does the total cross section. Certainly, adding a constant 
to V(r) does not affect the size of the cross section at the peak. Thus, we Fourier transformed 
Vjkpt numerically from momentum space to coordinate space and plot the derivatives of the 
potentials in Fig. @(a). To demonstrate the difference of the attractive forces, we show the 
difference of the derivatives of the two potentials, 

AF(r) = dVjKFT ^ _ ^SFHMNtr) 
dr dr 

(solid line) in Fig. 0(b). We confirm that AF(r) > holds in the region probed by the 
toponium states, r ~ 0.03-0.1 GeV -1 . One also sees that both potentials have a common slope 
at r > 0.4 GeV -1 because of the severe constraints from the bottomonium and charmonium 
data. The kink seen in the figure is due to a discontinuity of g^Vsfhmn/g^" 2 located at the 
continuation point, r = r c . 

In order to compare the asymptotic behavior of the potentials more clearly, we plot in 
Fig. [|(a) the coordinate-space effective couplings defined by 

a v (l/r) = (-C J r/r)- 1 y(r) (3) 

for Vjkpt (r) and Vsfhmn(^) as solid and dashed lines, respectively. Contrary to our expectation, 
the difference of the couplings exceeds 3% even at very short distances, 1/r ~ 100 GeV. 

Naturally the question arises: Why is there such a large discrepancy between the potential 
constructed in momentum space and that constructed in coordinate space? To answer this 
question, let us examine a relation connecting the effective coupling in coordinate space, defined 
by Eq. (|3|), and the effective coupling in momentum space, defined from the momentum-space 
potential as 

a y (q) = (-4nC F /q 2 y 1 V( q ). (4) 

The relation is derived from the renormalization group equation of ay(q) and exact to all 
orders. In the asymptotic region where the couplings are small, it can be given in the form of 
an asymptotic series JI3| , which reads numerically 

av(l/r) = «v + 1-225 a% + 5.596 ay + 32.202 oc% + . . . (5) 



' In this paper we are not concerned with those differences of the cross sections which can be absorbed into 
an additive constant to the potential, or equivalently, into a redefinition of the top quark mass. 
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Figure 2: (a) Comparison of the derivatives of the potentials vs. r for (M§) = 0.115: 
dVjKPT /dr (solid line) and o^Vsfhmn /dr (dashed line), (b) Difference of the derivative of the 
potentials vs. r. The solid line shows AF(r) = dV^w/dr — g?Vsfhmn / dr , and the dotted line 
shows AF(r) = dVj K p>r/dr — dV ncw /dr. 
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Figure 3: (a) Comparison of the coordinate-space effective charges defined from Vjkpt (solid) 
and from Vsfhmn (dashed line), (b) Comparison of the coordinate-space effective charges 
defined from the various terms of Eq. (Eh. See the text for the description of each curve. 
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for rif = 5. On the right-hand-side, ay — oty(q — e~ 7B /r). All terms which are written 
explicitly are determined from the known coefficients of the f3 function, 0q , (5\ , and 0% . At 
present, we can use the above relation consistently only at 0(a\) because we know the effective 
couplings only up to the next-to-next-to-leading order corrections in perturbative QCD, i.e. we 
know the relation between ay and a^ only up to O(a^). Due to this limitation, essentially, 
the effective coupling a v defined from Vsfhmn is the right-hand-side of the above equation 
truncated at the 0(a\) term, while ay defined from Vjkpt is the right-hand-side including 
all terms. Numerically, the 0(a,y) term and the 0(a\) term contribute as +1.5% and +1.2% 
corrections, respectively, for ay — 0.14 (corresponding to 1/r = 100 GeV). Therefore, these 
higher order terms indeed explain the difference of the effective couplings at small r. Fig. 0(b) 
shows several curves derived from the above relation: 

1. The solid line is ay (1/r) defined from Vjkpt- 

2. The dashed curve is ay + 1.225 a v , where ay = ay(q = e~ lE /r) is calculated using 
the perturbative prediction in momentum space. This curve is essentially the same as 
ay (1/r) defined from Vsfhmn ( r ), since it is the next-to-next-to- leading order perturbative 
prediction for the coordinate-space coupling at short distances. 

3. The dotted curve includes the next correction, 5.596 ay, which is in fact even larger than 
the 0(a\) term below 1/r ~ 30 GeV. 

4. The dash-dotted curve includes the 0(a\) term. 

We observe that the agreement of both sides of Eq. (|) becomes better as we include more 
terms at small r, while it becomes worse at large r on account of the asymptoticness of the 
series. From the purely perturbative point of view, the discrepancy between our two potentials, 
Vjkpt an d Vsfhmn, in the asymptotic region thus seems real, an indication of large higher 
order corrections.^ If the third-order correction to the potential will ever be computed in the 
future, the 0(a v ) term will be treated consistently and the difference will reduce by 1.5% at 
1/r ~ 100 GeV. 

The above 3% uncertainty of dy(l/r) at 1/r ~ 100 GeV provides a certain criterion for 
the present theoretical uncertainty of the ti cross section. In fact, it would already limit the 
theoretical accuracy of ay (1/r) at longer distances to be not better than 3%. If we combine 
this with a naive estimate <T pe ak oc oty, we expect a theoretical uncertainty of the peak cross 
section to be 5cr peak /crpeak ~ 6%. Therefore, the large discrepancy of the cross sections which 
we have seen turns out to be quite consistent with this estimated uncertainty. 

One is tempted to include one more term of the above series fl5|) to define a (new) coordinate- 
space potential despite our ignorance of the corresponding terms in the relation between ay 
and a^s, since this would apparently reduce the difference between the two effective couplings. 
In fact we did this exercise, but (to our surprise) it did not bring the cross section closer to the 
one calculated from the momentum-space potential Vjkpt- This cross section calculated from 

If we apply the same method (the relation between ay and ay) to estimate the size of the already known 
O(o.g) correction, we obtain 7t 2 /3q/3 = 193.4, which turns out to be a slight under- estimate of the true correction 
a 2 = 333.5 g| (n f = 5). 
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the potential V now (r), which incorporates the 0(a v ) term of Eq. (|5|), is shown as a dotted curve 
in Fig. 0(b). 

We may understand the reason why the cross section did not approach that of Vjkpt if we 
look at the difference of the "forces", AF(r) = d Vjkpt / dr — dV new / dr , shown as a dotted line in 
Fig. 0(b): it can be seen that, upon inclusion of the 0(a v ) term, the difference AF(r) decreased 
at small distances, r < 0.05 GeV -1 , as expected, whereas AF(r) increased at distances r > 
0.05 GeV -1 which is still in the range probed by the toponium states. It is due to a compensation 
between the decrease and increase of AF(r) that the normalization of the cross section scarcely 
changed. The increase of AF(r) at large distances results from the bad convergence of the 
asymptotic series, Eq. (||), for a large coupling, as we have already seen in Fig. |3](b). This fact 
indicates that we are no longer able to improve the agreement of the cross sections by including 
even higher order terms, as we are confronting the problem of asymptoticness of the series. 

Some indications can be obtained by looking into the nature of the perturbative expansion 
of each potential. Within our present knowledge of the static QCD potential, the perturbative 
series looks more convergent for the momentum-space potential than for the coordinate-space 
potential. To see this, one may compare the f3 functions of the effective couplings (the V-scheme 
couplings) in both spaces flO| . Numerically, the first three terms in the perturbative expansion 
read 

• (momentum-space coupling) 

^ 2 dav = _ Q 61Q1 q 2 _ Q 2UQ a s _ i.i98 a 4 + ( 6 ) 

(JjLL 

• (coordinate-space coupling) 

^ 2 dav = _ Q 61Q1 -2 _ Q2U9 -3 _ L945 (7) 

(JjUj 

for rif = 5. The first two coefficients are universal. The third coefficient depends on the scheme 
(the definition) of the coupling. As the third coefficients for the V-scheme couplings are quite 
large, the third term of the j3 function is comparable to the second term already for ay = 0.20 
and for ay = 0.13, respect ively.[] The difference of the third coefficients between momentum 
space and coordinate space originates from the 7t 2 /3q/3 term in Eq. (0), which comes from the 
Fourier transformation. (Compare Eqs. ( |T0| ) and (0).) Although the magnitude of the third 
coefficients is of the same order, in practice it makes a certain difference whether an apparent 
convergence is lost at ay = 0.20 or ay = 0.13 because there is a large scale difference between 
the two values. This indicates a worse convergence in coordinate space than in momentum 
space. 

§ The shift of the peak position to lower energy is caused mostly by a decrease of the constant Co in Eq. ( |l3| ) 
and not due to an increase of the attractive force. Since the effective coupling ay(l/r) runs faster for V new , the 
perturbative potential is connected to the intermediate-distance phenomenological potential at a deeper point. 

» This is the reason why we evolve the MS coupling instead of evolving the V-scheme couplings using their 
own P functions. Otherwise we would have lost the reasoning to keep the third term of the p function at a 
fairly large momentum/short distance. For comparison, the (3 function of the MS coupling for the same n/ is 
given by 

^ d*gg = Q 61Q1 a2 _ _ Q 2449 af _ _ Q9116 a 4_ + .... 
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Figure 4: Comparison of the cross sections for different choices of the scale. The upper 
three curves are for the momentum-space approach: \i = q (solid), \i = \/2q (dotted), and 
fj, = q/y/2 (dashed line). The lower three curve are for the coordinate-space approach: /i = 
/i 2 = exp(— 7b)/V (solid), /x = V2^ 2 (dotted), and /i = /i 2 / 'y/2 (dashed line). 

If we evolve the coordinate-space coupling ay using its own (3 function up to the third term, 
the coupling exhibits an infrared pole atl/r = A~2 GeV, which is an order of magnitude 
larger than Aj^g of the MS coupling. The asymptoticness of the series in Eq. ([5]) is closely 
related to the existence of this pole. In fact, one may estimate the uncertainty caused by the 
asymptoticness of the expansion to be 5ay(l/r) ~ Ar + (Ar) 2 + . . . JOJ. If we translate this to 
the uncertainty in the slope of the coordinate-space potential, we obtain SF(r) ~ A 2 . This is 
in good agreement with the discrepancy AF(r) ~ 1-4 GeV 2 in the region r > 0.05 GeV -1 (see 
Fig. |2|), where the usability of the asymptotic expansion is already limited to the first two or 
three terms. (For r < 0.05 GeV -1 , one may reduce the difference by including more terms.) 

It is interesting to examine the level of uncertainties within each of the momentum-space 
approach and the coordinate-space approach by itself. Fig. [| shows how the cross section 
changes when we vary the scale by a factor of 2 in each approach: from /i = qj a/2 to ji — \[2q 
in Eq. (28) of JTOf in the momentum-space approach (upper three curves), and from [i = /x 2 / v^2 
to fj, — v / 2/^2 m Eq. (44) of [|l(| in the coordinate-space approach (lower three curves). For 
the momentum-space approach, the variation of the cross section is about 2% at the peak and 
around 0.6% for cm. energies above threshold. Meanwhile in the coordinate-space approach, 
the variation of the cross section amounts to 0.9% at the peak and 0.8% at larger cm. energies. 
These results may be regarded as an internal consistency check for each approach and even as a 
sign for the stability of the theoretical predictions. Nevertheless one should keep in mind that 
the internal consistency is not the same as the accuracy of the theoretical predictions. Since 
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the 0{a 2 s ) corrections to the potential resulted in an unexpectedly large modification of the 
total cross section (see Figs. [| and |6]), it would be legitimate to consider each of our results as 
accurate only if the same method could estimate the size of the next-to-next-to-leading order 
correction reasonably well and hence if the cross section became considerably less sensitive to 
the scale variation after including this correction. This is not the case in our problem, however. 
The very existence of a large constant at the next-to-next-to leading order (0,2 in [ 10 1 , see 
also Eqs. fllOD and (|14])), which generates these large modifications, may indicate also large 
corrections at even higher orders. Furthermore, if we consider the full set of the fixed-order 
0(a 2 s ) corrections to the cross section near threshold they are larger in size and even more 
scale dependent than the corrections to the potential alone. The theoretical uncertainty may 
therefore be larger than indicated by the study in this paper. 

Still there may be some possibilities to reduce the difference between the momentum-space 
potential and the coordinate-space potential in the region probed by the toponium states. An 
obvious point to be improved is to remove the discontinuity of Vgp HMN (r) at r = r c . If we 
employed a smoother interpolation of the perturbative potential to the intermediate-distance 
potential, we would have a better agreement of the two cross sections, see Fig. 0(b). This 
tendency is expected due to the specific interpolation method adopted for VsfhmnM-Q However, 
one has to be careful with this argument. The slope of the potential in the intermediate-distance 
region is fixed by experimental data, which correspond to one fixed value of «j^g(M|) (= the 
true value in nature). We are interpolating the prediction of perturbative QCD, which obviously 
depends on our input value of ajjg(M|), to a phenomenological potential, which is independent 
of it. This means, we do expect a non-smooth transition for any value of a^g(M|) different 
from the true value. Moreover, if we want to extract the value of a^(M^) by comparing the 
theoretical predictions to the experimentally measured cross section, the sensitivity to a^g(M|) 
decreases if the predictions depend on the way we perform the interpolation. Ideally we would 
want to have an intermediate-distance potential as the prediction of QCD — necessarily non- 
perturbative — for a given input value of o^jg(M|). 

It would be important to understand the problem of the difference in the potentials also 
in momentum space, at least as much as we do in coordinate space presently. We have not 
done this analysis so far because of the difficulty in the numerical Fourier transformation of the 
potential from coordinate space to momentum space. 



Summary 

• There is a difference between the potential constructed in momentum space and that 
constructed in coordinate space even at a fairly short- distance, 1/r ~ 100 GeV. The 
difference can be understood within the framework of perturbative QCD. We already 
know that there is a large correction at 0(af) in the relation between the two potentials, 
although a consistent treatment is not possible until the full 0(af) corrections to the 
QCD potential are calculated. 

" Vsfhmn(?") matches the perturbative potential exactly up to a vicinity of the infrared pole. Since the rapid 
acceleration of the running of ay(X/r) towards the pole tends to amplify the deviation from Vjkpt(^), if we 
employed an analytic regularization of the pole while keeping the potential to approximate the perturbative 
potential at short distances, the running of dy(l/r) should be tamed, and hence the potential should come 
closer to Vjkpt(^). 
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• The above difference at short distances provides a criterion for our present theoretical 
accuracy of the tt cross section, &x pe ak/0"pcak > 6%. 

• In addition, it seems that we are confronting the problem of the asymptoticness of the 
perturbative series in the calculation of the ti cross section, as the top quarks do not 
probe a region which is sufficiently deep in the potential. We may not be able to improve 
our theoretical precision even if the higher order corrections are calculated in perturbative 
QCD. 

• We may, however, discuss which of the two approaches gives a more favorable result 
theoretically Up to the second-order corrections, the perturbative series looks better 
convergent for the momentum-space potential than for the coordinate-space potential. 
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A The momentum-space potential 

It seems appropriate to describe the potential Vjkpt used in the present analysis in some more 
detail. This potential is very similar to the potential Vjkt described in and used in all later 
numerical studies within the momentum-space framework. It includes, however, the next-to- 



next-to-leading order terms from ||10|| . The momentum-space potential can be written as 

,2 



VW(q) = K)(qcut) ■ (2vr) 3 5(q) - 4ttC f ajK ^ q I , ( 8 ) 



The effective coupling ckjkpt is defined to coincide with the two-loop perturbative prediction 
for large momenta, to be Richardson- like for small momenta, and to simply interpolate between 
these two shapes in some intermediate range: 

!av, P crt(q 2 ), |q| > qi = 5 GeV 

«Rich(q 2 ), |q| < q2 = 2 GeV (9) 

aRich(q 2 ) + (av, P crt(q 2 ) - aRich(q?)), q2 < |q| < qi- 

The intermediate regime is only introduced to obtain a smoother transition between the small 
and large momentum parts, respectively. 

The first difference between the updated potential and the former version is the fact that 
we are now able to use the full two- loop expression for the perturbative part, 

,2\ _ „ t„2\ ( 1 i „ a Ms"( c l 2 ) _j_ a ^Msi^ >X 

47T 2 V 47T 



«V, P e rt (q 2 ) = «Ms(q 2 ) l + ai ^^ + a 2 p^) , (10) 



with the coefficients ai and a 2 given in |TIJ. As the 6-quark threshold is neglected, rif = 5 is 
set throughout in the evolution of the MS-coupling, which can now consistently be performed 
at three-loop accuracy. 
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A Richardson-like behavior for small momenta is chosen since the Richardson potential |14 
is known to describe the charmonium and bottomonium spectra fairly well. A pure Richardson 
form, however, would lead to severe numerical problems. Hence the ansatz has to be modified 
slightly by introducing two "subtraction terms" , 

. 9 . An f 1 Al Al \ . . 

aRich q 2 = -5-7 -r — -yr - -f + 2 _i_ 2 11 

with Ar = 400 MeV. The first subtraction regulates the divergent behavior for |q| — > 0, 
the second subtraction is designed to reduce the modification introduced through the first to a 
minimum. Without the second additional term, the linear part of the position-space Richardson 
potential would be removed completely, whereas with it the first subtraction is cancelled for 
q 2 ^ Qcut > an d thus a big part of the confining potential is kept. It thus seems desirable to 
choose the parameter q cut small, but evidently it cannot be put to zero to really recover the 
pure Richardson potential. However, the linear part of the potential plays practically no role 
for the tt-system as will be demonstrated below. The exact value of q cut is therefore relatively 
unimportant and the adopted value q cut = 50 MeV results in both numerical efficiency and 
speed of the program and a fairly good accuracy of the predictions. 

The constant Vo(q cu t) in Eq. @ is to some extent an arbitrary parameter. Different choices 
of Vo(q cu t) reflect the ambiguity in the definition of the pole masses for confined quarks. For 
VjKPT(q) the choice 

47rC F Al , . 

Vb(qcut) = 777 12 

#o(n f = 3) q cut 

is used. It leads to a Richardson-like potential that depends only weakly on the parameter q cut 
and coincides with the true Richardson potential in position space in the limit q cu t — > 0. With 
this potential one obtains for the pole mass of the b quark m& = 4.84 GeV. The choice of Vq is 
the second difference to the potential used in earlier works, where the constant Vq was fixed by 
the condition V JKT (r = 1 GeV -1 ) = -1/4 GeV leading to m b = 4.7 GeV. 

Fig. |5] shows a comparison of the total cross section for tt-production near threshold as a 
function of energy using the "old" potential, i.e. the one as described but using the one-loop 
formula for expert ( an d correspondingly the two-loop evolution for the MS-coupling) and fixing 
Vo through V(r — 1 GeV^ 1 ) = — 1/4 GeV (dash-dotted line), and the "new" potential (solid 
line). There are two changes: first, the inclusion of the two-loop correction to the perturbative 
potential increases the strength of the attractive interaction between t and t, and thus leads 
to an increase in the cross section. This is nicely demonstrated by the dashed curve, which 
corresponds to the inclusion of the two- loop potential, but the old choice of Vq. Second, the 
modified choice for Vq leads to a small shift of about 300 MeV in the energy scale, which is 
just the difference between the two Vq. The dotted curve has been included to once again 
demonstrate that the tt system is quite insensitve to the long-range part of the potential: this 
curve corresponds to the choice q cut = oo, i.e. to completely removing the linear part from the 
potential and setting Vq = 0. 

B The coordinate-space potential 

The short- distance part of the coordinate-space potential is given by the next-to-next-to-leading 



order static QCD potential in position space UlOl , whereas its form in the intermediate- and 
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Figure 5: Comparison of the total cross section (normalized to R) for tt-production as a func- 
tion of E = ^/s — 2m t for m t = 175 GeV, = 0.115, and T t = 1.427 GeV, using the different 
potentials described in the text. The dash-dotted line corresponds to the model of older anal- 
yses, the dashed line to the inclusion of the two- loop correction to the perturbative part. The 
solid line shows the prediction of the new model, which differs from the dashed line only by a 
different definition of the energy. The dotted line corresponds to the new potential with the 
linear part in position space removed. 



long-distance region is determined phenomenologically. We thus have 



SFHMN 



J") 



v; 



pert I, 



at r < r, 



c + C\ log(r/r ) exp(— r/r\) + ar at r > r, 



C ) 

c ■ 



(13) 



Here. 



^pert(r) = -Ci 



1 + Oi 



Air 



+ (a 2 + 



An 



(14) 



represents the coordinate-space potential in the second scheme, = exp(— r fE)/r. The coeffi- 
cients a\ and a2 are the same as in the momentum-space potential, and (3q = (IICa — 4Tprif)/3. 
See Ref. for details.Q 

The values of the phenomeno logical parameters r , r±, a and C\ are taken from Ref. |3J and 



** We evolve the MS-coupling ct^iji) by solving the three-loop renormalization group equation numerically 
for a given initial value at /i = Mz, whereas an approximate solution to the renormalization group equation is 
used in |l(J. 
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are tuned to reproduce bottomonium and charmonium data well: 

r = 0.2350 GeV" 1 

ri = 3.745 GeV" 1 ( , 

a = 0.3565 GeV~ 2 ^ ' 

ci = 0.8789 GeV 

We fix c and r c by requiring that both the potential Vsfhmn (?") and its first derivative are 
continuous at r = r c . For example, r c = 0.2526 GeV -1 and cq = —1.972 GeV for a^g(Mf ) = 
0.115. 

This potential is an improved version of the potential proposed in ]3[ by including the next- 
to-next-to leading order terms to the short- distance QCD potential. We compare the cross 
sections calculated from the present version and from the old version in Fig. |G] for a^jg(M§) = 
0.115. 



R(E) 




4 

E [GeV] 



Figure 6: Comparison of the total cross section (normalized to R) as a function of E = ^—2mt 
for m t = 175 GeV, = 0.115 and and T t = 1.427 GeV using the old and the present versions 
of the coordinate-space potential Vsfhmn- The dash-dotted line corresponds to the potential 
described in ||. The solid line shows the prediction of the present potential, Eq. (p~3|). 



13 



References 

[1] M.J. Strassler and M.E. Peskin, Phys. Rev. D43 (1991) 1500. 

[2] V.S. Fadin and V.A. Khoze, JETP Lett. 46 (1987) 525; Sov. J. Nucl. Phys. 48 (1988) 309. 

[3] Y. Sumino, K. Fujii, K. Hagiwara, H. Murayama and C.-K. Ng, Phys. Rev. D47 (1993) 
56. 

[4] M. Jezabek, J.H. Kiihn and T. Teubner, Z. Phys. C56 (1992) 653. 
[5] M. Jezabek and T. Teubner, Z. Phys. C59 (1993) 669. 



[6] M. Peter and Y. Sumino, University of Karlsruhe Report No. TTP97-27, [hep-ph/9708223 
to appear in Phys. Rev. D. 

[7] A.H. Hoang, Phys. Rev. D56 (1997) 5851. 

A.H. Hoang, P. Labelle and S.M. Zebarjad, Phys. Rev. Lett. 79 (1997) 3387. 

[9] G. Adkins, R.N. Fall and P.M. Mitrikov, Phys. Rev. Lett. 79 (1997) 3383. 

[10] M. Peter, Phys. Rev. Lett. 78 (1997) 602; Nucl. Phys. B501 (1997) 471. 

[11] A. Czarnecki and K. Melnikov, University of Karlsruhe Report No. TTP97-54, hep- 



| ph/9712222| , and references therein. 

[12] A. Hoang and T. Teubner, University of California, San Diego Report No. UCSD/PTH 
98-01, DESY 98-008, |hep-ph/9801397| . 

[13] M. Jezabek, M. Peter and Y. Sumino, in preparation. 

[14] J.L. Richardson, Phys. Lett. B82 (1979) 272. 



14 



